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Abstract. Previously a many-body coherent potential approximation (CPA) was 
used to study the double exchange (DE) model with quantum local spins S, both for 
S = 1/2 and for general S in the paramagnetic state. This approximation, exact 
in the atomic limit, was considered to be a many-body extension of Kubo's one- 
electron dynamical CPA for the DE model. We now extend our CPA treatment to 
the case of general S and spin polarization. We show that Kubo's one-electron CPA is 
always recovered in the empty-band limit and that our CPA is equivalent to dynamical 
mean field theory in the classical spin limit. We then solve our CPA equations self- 
consistently to obtain the static magnetic susceptibility x in the strong-coupling limit. 
As in the case of the CPA for the Hubbard model we find unphysical behaviour in x 
at half-filling and no magnetic transition for any finite S. We identify the reason for 
this failure of our approximation and propose a modification which gives the correct 
Curie-law behaviour of x at half-filling and a transition to ferromagnetism for all S. 



PACS numbers: 75.20.Hr, 75.30.Kz, 75.40. Gb 



1. Introduction 

Recently there has been much interest in the perovskite manganite compounds 
Ti_a.D3-Mn03 where T and D are trivalent and divalent cations respectively. 
These exhibit a rich variety of phases including charge, orbital, ferromagnetic and 
antiferromagnetic ordering |]. Of particular interest is Lai_2;Caa;Mn03 with x ~ 0.3; 
in this compound ferromagnetic-paramagnetic and metal-insulator transitions occur 
together, and for temperatures near the critical temperature an applied magnetic field 
causes a very large reduction in electrical resistance: this is the phenomenon known as 
colossal magnetoresistance (OMR). 

The physically relevant electrons in the manganites are those from the Mn 3d levels, 
which are split by the approximately cubic crystal field into triply degenerate t2g levels 
and higher energy doubly degenerate Cg levels. Occupied Cg levels are further split into 
two nondegenerate levels by the Jahn- Teller effect. Electrons from the eg levels are able 
to hop between Mn sites via the O atoms, forming a narrow conduction band, but those 
from the t2g levels are localised. There is a strong Hund's rule coupling on the Mn 
sites, so the t2g electrons are usually modelled as local S = 3/2 spins ferromagnetically 
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coupled to the itinerant Cg electrons. At a; = there is one Cg electron per site so the 
system is a Mott insulator. Doping by D atoms produces holes in the Cg band which 
enables conduction to occur. 

The simplest model for the CMR compounds, which neglects the Cg degeneracy and 
any coupling to phonon modes, is Zener's double exchange (DE) model 

h=y: ua.c,. - j e ■ - e (1) 

ija i i 

where i and j are Mn sites, c^^ (c|^) is a cr-spin conduction electron annihilation 
(creation) operator, Sj is a local spin operator, cTj is a conduction electron spin operator, 
Lf = Sf + erf is the 2;-component of the total angular momentum on a site, Uj is the 
hopping integral with discrete Fourier transform tk, > is the Hund's rule coupling 
constant, and h = g^sB is the Zeeman coupling strength, B being the applied magnetic 
field. The number of conduction electrons per atom n is assumed to be given by 
n = 1—x. The idea of the DE model is that hopping of Cg electrons between neighbouring 
sites is easier if the local spins on the sites are parallel, so an effective ferromagnetic 
coupling between the local spins is induced by the conduction electrons lowering their 
kinetic energy. Double exchange coupling differs from conventional Heisenberg coupling 
by being (for classical local spins) of the form cos{9ij/2) rather than cos(^ij), where 9ij 
is the angle between the i- and j-site local spins. 

According to Millis et al ||, g the CMR effect arises from a competition 
between double exchange coupling, which produces a tendency towards the conducting 
ferromagnetic state, and strong coupling to phonon modes, which tends to localise 
the electrons via self-trapping. In a previous paper we confirmed that the DE 
model above cannot account for the very high resistivity of the paramagnetic state. 
Moreover experiments show that coupling to the crystal lattice is important P, 0, §]. 
In this paper however we will complete our study of the simple DE model, aiming to 
understand the purely electronic properties of CMR systems modelled by (|I]) before 
tackling more realistic and complicated models. We concentrate particularly on the 
magnetic properties. 

In 1^ we derived an approximation for the one-electron Green function which was 



based on Hubbard's scattering correction approximation for the Hubbard model [10 
In the Hubbard model this approximation is derived by decoupling the Green function 
equations of motion according to an alloy analogy in which electrons of one spin are 
frozen whilst the Green function for those of the opposite spin is calculated. For 
the finite S DE model this approach is complicated by the possibility of dynamic 
spin scattering — conduction electrons exchanging angular momentum with the local 
electrons. We obtained an approximation which, like Hubbard's, was exact in the 
atomic limit for all band filling. Since Hubbard's approximation is equivalent to the 
coherent potential approximation (CPA) in the alloy analogy for the Hubbard model, 
and our approximation reduces to a one-electron dynamical CPA due to Kubo [|ll[ 



m 



the empty-band limit of the DE model we regard our approximation as a many-body 
extension of the CPA. 

In P] we concentrated on simple cases in which the CPA system of equations 
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of motion closed easily, and calculated the electronic structure and resistivity of the 
paramagnetic state. In section Q we formulate and solve the general CPA equations. 
In section |^ we then calculate the static magnetic susceptibility self-consistently within 
the CPA. In section ^ we compare our CPA in the classical spin limit with dynamical 
mean field theory, and in section |^ the CPA is modified so as to improve the behaviour 
of the susceptibility. A summary and outlook are given in section |[ 



2. Solution of the CPA equations 

In this section we will use equation of motion (EOM) decoupling approximations to 
derive an expression for the one-electron Green function = ((Cj^;c]^)) of the 
DE model. In a previous paper [^j we obtained in the special cases of zero-field 
paramagnetism and saturated ferromagnetism for all values of S, but only considered 
the case of arbitrary magnetisation for S = 1/2. Here we extend this previous treatment 
to the case of general 5* and spin polarization. The decoupling approximations used are 
direct extensions of those used in P], which were in turn generalisations of Hubbard's 



scattering correction approximation ||10| for the Hubbard model. 



We split the Green functions into components G^J = G^J^+G^J where G^" describes 
propagation via singly (a = — ) and doubly (a = +) occupied sites: G^^" = {{nfc^^ ; c]^.)) 
where n+ = n-, = 1 - n-, n+ = n-^, and ra" = 1 - n^^. Here n^^ = c^Ci^ 
and rii = rii^ + n^j^. Now the cases considered in were chosen so that the system 
of EOM could be closed using only the Green functions GJ^"", {{Sfnfc^^]c^j^)), and 
{{Sj^'^n°'Ci_^ ] Cj^)) where S~'^ = 5*^" and for cr =t, i respectively; here we will also 
need the Green functions (((Sf )™?T,°Cj^ ; cjg.)) and (((5'j^)™~^S'~°"ra"Cj_^ ; c]^)) for m > 1 
in order to close the system. It is simplest to work with the Green functions 

SrW = {{e''^Kc,^;cl)) (2a) 
= {{e''^ Srntc^^^ ; 4)) , {2b) 

which are in generating function form so that d"^S^J°'/dX^\x=o = (((S"/ )™n"Cj^ ; c]^.)) 
and a— iT^^V^A— 1|a=o = {{iSn"'-'S-nfc,^ ; 4)). 

For notational simplicity we will work for a =1; the cr =| EOM can be obtained 
by making the replacements Qo- i— > Ci-a, c\„ ^ cl_^, S- t— > —Sf, Sf i— > S^, 
h \—>- —h and A t-^ —A. Recall that with the fermionic definition of Green functions, 
{{A;C)), = -tj^dt exp{iet){{A{t),C}), the EOM is 

e{{A;C))^={{A,C}) + {{[A,H];C))^. (3) 

We use the fact that S^J°'{X) and T^-'"(A) are in generating function form to write their 
exact EOM in the form 



^ ^ (A,e) + i 



k 
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for 5'|-'"(A, e) and 



+ 7^ 1 1 ^r(A, e) + -e-^^- {SiS + l) + a-- — ] SriX, e) 




k 

+ {{^''^S-[n-,H,]c,,;cl,))^ (5) 

for T|-'°'(A,e). Here Hq = Y^ijaUjcl^Cj^ is the kinetic part of the Hamiltonian and 
we have dropped the site indices in the expectations, assuming the system to be in a 
homogeneous phase (i.e. we will not consider antiferromagnetism) . If one works directly 
with EOM for the Green functions )"^<Ci^ ; c]^)) and (((5f )"^i5r<Cij, ; c]^)), 
as in one is faced with 45* + 1 algebraic equations for general S. Here, using 
the generating functions, these algebraic equations are reduced to just two differential 
equations where differentiation with respect to A corresponds to the coupling between 
the different algebraic equations. 

These EOM form a closed system apart from the undetermined Green functions 
on the right-hand sides which correspond to the effects of hopping. The decoupling 
procedure entails making approximations for these kinetic terms which close the system 
of equations; since these terms are proportional to t this procedure is exact in the 
atomic limit tij -^0. As mentioned in the introduction the idea of the alloy analogy 
is to neglect the effects of the kinetic part Hq of the Hamiltonian on electrons of one 
spin whilst considering the propagation of an electron of the opposite spin. Accordingly 
we neglect the final terms of (^) and (^) since the occupation number operators in 
these Green functions are considered to be frozen in time. It remains to find closed 
approximations for Ea; ^ifc((exp(AS'f X^c^p c]|))e and Y.kUk{{exY>{XSI)S~nf^c,,^]c]^))^. 
The derivation of scattering correction approximations for these terms is identical to 
that of apart from the occurrence of the factor exp(AS'f) and will not be repeated 
here. In fact we use 




(6a) 

JT(e)Sf (A,e) (66) 



+ J^[e + h){{e''tS-r^^c,^■c]^) + J^{e) UaU'^' S'aA c,^;c]^) (6c) 



—a 




t,uG\\e) - Me)G]\e) ]+J^{e + /i)Tf "(A, e){Qd) 



k 
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where Jo-(e) = e — So.(e) — G(y{t) . Here So-(e) is the self-energy, local within this 
approximation, and Ga(e) = G^*(e) is the local component of the Green function. Jcr(e) 
contains the effects of coherent propagation of the electron as a a-spin from site i back 
to site i via paths avoiding the site at intermediate stages |T^. Note that if an |-spin 



electron of energy e becomes a |-spin by exchanging angular momentum with a local 
spin it must then propagate at energy e + /i, hence the occurrence of Jj(e + h) in (|6^ ) 
and ( |6^ ) above. It may be seen that approximations and close the system of 
equations and (^. No further approximations are made. 

For convenience we now define -Eo-(e) = ^ ~ Jai^) , which will later be related to the 
Weiss function of dynamical mean field theory [Q, E^{e) = E„{e + hS^i) +ah/2, which 
puts the energy shift effects of the magnetic field into E^-, and 

= S., +Y^U,G'j{e) - Me)G^^{e). (7) 



We make the above approximations so that and become 
J/2e-''-^[S{S + l) + a-§-,-£] Ef{e) - J/2 + £ 



This is a coupled pair of linear differential equations (with respect to A) of first and 
second order respectively. The first equation is used to eliminate T|"'°(A,e) in terms of 

Tf "(A, e) = |e-^^- (^Af (e) (e^^^nj) - (e)5f (A, e) - ^-^S^iK e)) • (9) 

Substituting into the second equation we obtain a first order equation for S'|-'"(A, e), 
d ( J/2SiS+l)-E^ie)i2/JE^ie) + a) \ _ 

/(2/JEf(e) + a-^^)e^^'n5* + ae^(^"++^')S-a+\ .. 

This equation is of the form 

^ + PS{X) = R{X), (11) 

where P is independent of A, which has the solution 

^(A) =Ce-^^ + e-^^y'^dAe^^i?(A) (12) 

where C is a constant of integration. Furthermore by inserting local spin projection 
operators into definition (p^) of S'|-'"(A,e) it may be seen that S'|-'"(A,e) = 
J2m=-s '^m^i^) ^^Pi''^X) where a^"(e) is independent of A. Since in general P is not 
an integer we must have C = 0. Hence we find 



Sr{X,e) 



E\{e) + J/2(a - S'))e^^^nl - {J/2 6a-)e^^^ 



(Ef (e) + J5V2) (^f (e) + J/2{a - S^)) - jy4{S + aS^){S + 1 - aS^ 
(J/2 5,+)e^(^^+i)5-a+ 



(Ef (e) + J/2{S' + 1)) (Ef (e) - JS'/2) - J^/A{S - S'){S + I + 



Xne) (13) 



where we adopt the convention that a quotient of operators, B divided by A, means 
A^^B. Substituting into (||) we also obtain 

..^ /e^(^'-^-+)(j/2(^^ + aS){S' - a{S + + (e) + JS'/2)S-a+) 

^ \ (^f (e) + JS^/2) (E\{e) + J/2(a - 5^)) - JV4(5 + aS'){S + 1 - a5^) 

-7 \, ^- ^-^ )Ay(e). (14) 

(Ef (e) + J/2(l + ^-)) (Ef (e) - - JV4(S - + 1 + 5^) / ^ 

It is easy to check that these expressions reduce to those obtained in ^ in the appropriate 
cases. 

We are mostly interested in the local components of the Green functions, S'^ = S*"" 
and = T|*", and in the following we drop site indices. Note that A^*(e) = 1; this 
follows from the relation A^ (e) = G^ (e)/Go-(e) which is easy to obtain from the definition 
(0) of A^(e) using Fourier transforms and the locality of T,„{e), as in section 3.1 of [0. 

Since Ej^ is a functional of G^, (|TB|) with i = j determines G^, in principle, in terms 
of the expectations (exp(AS'^)), (exp(AS'^)no-), and (exp(AS'^)S'^cr^) as functions of A. 
These expectations must be evaluated self-consistently in terms of 5*" (A) and T°(A). The 
equation for G^j simplifies if we approximate the physical (simple cubic tight-binding) 
bare density of states (DOS) by an elliptic DOS, ^(e) = 2VW^ - / (ttW) where W is 
the half-bandwidth; for this DOS it may be shown that E„{e) = e — W'^/4:Go-{e) so that 
Erj is an explicit function of Ga- We now define the functional 

IM = - /i) ^7(e) (15) 

where / is the Fermi function, fi is the chemical potential, and 7 is the anticlockwise 
contour lying just below and just above the real axis. This functional is useful owing to 
the usual sum rule 

1[{{A;G))] = {CA), (16) 

from which it follows (exactly) that 

I[S:{X)] = {e''\-^n^) (17a) 

/[T;(A)] = (e"^^5-V-)5„_. (176) 

The expectations (exp(AS'^)no-) and (exp(AS'^)S'^(j^) can therefore be obtained directly 
from the sum rule as (exp(AS'^)no.) = I[S^{X)], where S^{X) = J^a^aW^ 
{exp{XS^)S~a^) = /[T|(A)]. However (exp(AS'^)) must be obtained indirectly, in 
principle, by solving the system of equations obtained by applying the I-functional 
to (|13D and ([l^) . Knowledge of (exp(AS'^)) is equivalent to knowledge of P{S^), the 



probability distribution function for local spins, and it is not clear that this quantity 
will be accurately obtained from the (approximate) single-electron Green functions that 
we have considered here. In fact, as will be seen later, the self-consistent determination 
of this quantity causes problems with our CPA. 

We now specialise to three important cases of particular interest: n = 0, J = 00, 
and = 00, in which (|T^) and ([1^ simplify considerably. The results below generalise 
previous work H] which was restricted, for general magnetisation, to the case 5 = 1/2. 
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2.1. The empty band limit 



In |Tl| Kubo used a one-electron dynamical CPA to derive an expression for valid 
in the low-density limit n ^ 0. From ([T^) with A = we calculate in this limit as 



/ E'^ - J/2(l + S') \ , , 

^ ^^1^ + J/2 5^) (^f - J/2 (1 + 5-)) - JV4 (5 + 1 + 5^) (5 - 5^) / 

This is equivalent to Kubo's equation for G-^ so our decoupling approximation is indeed 
a many-body extension of the CPA. 



2.2. The strong coupling limit 

In the physical systems for which the double exchange model was introduced J ^ tij 
and < n < 1. In this situation the chemical potential lies in the lowest band near 
— JS/2, so we shift the energy origin, E^]; — JS/2, and let J ^ oo. Equations 

(p!3|) and (p!4[) then become 

{S+1 + S')nl + S-(x~ 




{S + 1 + S') E^ + {S-S' 



{S + 1 + S')n-^ + S-a+ \ 
' ^^\s + l + S^)E\^ + {S-S^)E\l ^^^^^ 

2.3. The classical spin limit 

In dynamical mean field theory (DMFT), a local approximation exact in the infinite 



dimensional limit 13 1, the double exchange model can be solved exactly in the classical 



spin limit S* — > oo in which the local spins can be treated as static [0]. Since our 
approximation is also local it is interesting to compare our results with DMFT in this 
limit. We let — > oo in (|TEp and (0), scaling J, A, h and T^(A) as l/S", and obtain 

{E^ - J/2 S')n^^ + a J/2 5"^+ \ 




(Et + ^/2 S^){E^ - J/2 S^) - JV4(1 - {S^Y) 

J/2{{S-'f - l)nj - a{E^ + J/2 S')S~a+ 
(Et + JI'^S^){E^ - J/2^^) - JV4(1 - {S^Y) 



(20 tt) 
(206) 



where by and S~ we mean S'^/S' and S^/S respectively. Note that in this limit 
Ej^ = E„. In section ^ below we will derive these Green functions within DMFT for 
comparison, and it will be found that our CPA agrees with DMFT for S*" and T^, but 
not for T^. 



3. Self-consistent CPA susceptibility 

In this section we calculate the static magnetic susceptibility x of the zero-field 
paramagnetic state. For simplicity we specialise to the strong coupling limit J = oo most 
favourable to ferromagnetism and use the elliptic bare DOS mentioned in the previous 
section. We drop the spin suffices on zero-field paramagnetic state quantities and define 



6A to be the first order deviation of any quantity A from its value in tliis state. Tlius 
6 A is proportional to the applied magnetic field B or equivalently to S{L^) = 6{S^ + cr^). 
We proceed by calculating S{L^) in terms of the Zeeman energy h = q^bB and using 
X = {gi^Byhmh-^o{5{L')/h). 

We first derive a couple of useful identities. Equations ( |19fl| ) and ( |19/] ) imply that 
To-(A) = {S — d/dX)Sa{X)- We apply the sum rule ([l6|)-( pT^ ) to this relation and, using 
the fact that all the self-consistently determined expectations are real, obtain the rather 
obvious results 

[e'''S-a) = S/2 (e'^'n) (21a) 
(S'n) = 2S {a') . {21b) 

These serve as a check on our approximation and will later be used to manipulate 
expectations. 

The Green functions are especially simple in the (J = oo) zero-field paramagnetic 
state: from ( |19q| ) with A = and a = — we have 

_ {S + l-n/2)/{2S + l) 
^^'^ - m ' ^ ^ 

which corresponds to a band of weight (5 + 1 — n/2)/{2S + 1) per spin. In the elliptic 
DOS case, for which E{e) = e — /4:G{e), we can easily solve this equation to obtain 

2 



G(6) = — (e-Ve^-W^M (23) 



where W = W y {S + 1 — n/2)/{2S + 1) is the half-width of the renormalised band, 
which is also elliptic. 

We expand the denominator of ( |19q| ) in powers of 5*^, 



r=o ((5 + 1)E!^ + SE\ 



Owing to the presence of the factor E^ — E^, which is zero in the /i = paramagnetic 
state, only the r = and r = 1 terms contribute to 5S^{\). From we find 

, , 5EHd/d\ -S)- 5EHdld\ + S + 1) , . , \ 

^^t(A) = ^^^^ ^Wt)!^ ^(^ ((5+l + 5>,- + 5-a+)) 



6{e^^ {(S +1 + S')n7 + S-a^ 

^ i2S,l)E ^ P^' 

where the paramagnetic state value, which may easily be evaluated, is used for the first 
expectation. Now it may be shown that 6fi = so that Sl[A] = l[SA] for any A. We 
calculate 5I[(9/9A(S'|(A) + S'j^(A))]a=o from (pSf), and using ( pioD and (|21£|) and the sum 
rule ( |17a| ) find 

6 {S^n^ +n^))=2S6 (a^) = ^ ^ ^_ {{S + 1)6 {U) - {23 + 1)6 (a^)) 



2S{S +l-n/2) 
3(25 + 1)2 



6E'l - 2{S + 1)6E'^^ 
^2 



(26) 
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Since E^{e) = e + h/2- W^AG^ie + h5,^) we have 

= h/2 - W^/A {5G„ + h5,^ G') (27) 



where G'{e) = dG(e)/de. By setting A = in (^) and using (^, (^) and (^ we 
obtain 

S E + h {W/Wr{SW'G'/6 -{S + 1/2)) 
" (25 + 1)^2- (2/3 5 + 1)W^74 ■ ^ 

Note that since 6G^ oc a spectral weight is transferred between the different spin bands 
at constant energy, so = as mentioned above. Equation ( pH]) illustrates how S{L^), 
and hence S{S^), is determined indirectly in terms of the Green functions rather than 
by a direct sum rule of the type ( P.7g| )-( [T7^ . 

We use (pTl) and ( pH]) to eliminate the 6E^s from ( p6D in terms of h and By 
applying I to (pSj) we obtain an expression for 5{a^). We then eliminate 5(cr^) between 
this equation and ( p6D and solve the resulting equation for x = (fl'/^B)^<^(-^^)/^, 



X = - 



\2 



2^ (4/3 5(5+l) + (4/3 5+l)n)i!^ (29a) 



g = i 



(296) 



where i/^ = (2/3 5* + l)/(25' + 1). Note that for S = 1/2 this expression for x 
does reduce to the one given (for the S = 1/2 case) in |jl5[. We can simplify 
Q by changing variables to ^ = W^G{e)/{2W) so that de = 2W /W^ G'{e)-^dz = 
W/2 (1 - {2W/WyG{e)-^)dz = W/2 (1 - z-^)dz, and from the functional form (0) 
of G it may be seen that the contour 7 for e becomes —7', the clockwise unit circle, for 
z. Hence 

where /i = n/W and /(e) = f{We). The same change of variables can be used to show 
that 

no'i^-Zl^fi^-,]. (31, 



Now a second-order transition to ferromagnetism corresponds to a divergence in 
X, i-e. Q = 1. From (|28| ) with /i = and (|29^ it may be seen, using the sum 



rule, that 6{L')Q/{2S + 1) = l[SG^] = 5{a') (for h = 0). The equation Q = 1 for 
a zero field magnetic transition is therefore equivalent to the consistency condition 
25(0-^) = 5(5"^ + cr^)/(S' + 1/2), which certainly holds at = 1 but not at n = 0. 
Integral (^) can be evaluated analytically in the limits of zero and infinite temperature 
where the Fermi function is of a simple form and the results are plotted in figure 0. It is 
clear that within the CPA there is no magnetic transition for < < 1, as is the case 
for the CPA for the Hubbard model [16]. This appears to be a considerable drawback 
of our approximation given that the ferromagnetic-paramagnetic phase transition is a 
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Figure 1. The Q function versus filling n at temperature T — and T ~ oo. 

major reason for interest in the double exchange model. We will propose a method for 
circumventing this problem in section ^. 

We now consider the behaviour of x at n = 0, 1. In these cases electron hopping 
does not occur and the system consists of a lattice of free local moments of magnitude 
S' = S and 5* + 1/2 respectively, so we expect x to take the Curie law form, 
Xc = {qI^bY I3S' {S' + l)/3 where /3 = {k-QT)~^ is the inverse temperature. We calculate 
X in these cases by expanding the integrals in (|30| ) and ( |3TD in powers of n or 1 — n. At 
n = we find x = Xc with S" = S" as expected, but at n = 1 we find x = Xc0 (with 
S' = S + l/2) where 



L, dzexpil3iz + z-^)/2) 



L, dz exp{(3{z + z-i)/2) {z - z-^)/{z^ - u^) 



(32) 



and P = W(3. Now — 1 as /5 — so the Curie law is obtained at high temperature, 
but as /5 — > oo we find 8S'^/(3(2S' + 1)(4S' + 3)) 7^ 1 so Curie law behaviour 
does not extend over the whole temperature range. Note that (f){(3 = 00) 7^ 1 even 
at S* = 00 where, as will be shown in the next section, our Green function equations 
reduce to DMFT. The reasons for this unphysical behaviour, and a way to avoid it, 
are discussed in section |[ In figure |^ below we compare x at n = 1 and S = 00 with 
Xc and the low temperature asymptote Xc/3. In figure || we plot x^^ at S* = 00 and 
n = 0.75, comparing it with the Curie law and DMFT values, the DMFT plot only 
being displayed for the paramagnetic phase. 



4. Comparison with dynamical mean field theory at S = cxd 



In this section we will obtain equations for the 5" and Green functions within 
dynamical mean field theory (DMFT) at = 00 and compare them with equations 
(M and (M of the CPA. The DMFT single-site effective action S of the DE model 
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Figure 2. The S = J = oo magnetic susceptibility x for n = I compared with its low 
temperature asymptote and the Curie law. 



IS 



a 



Jo Jq „ 2 Jo L 



+ 5+(r)a-(r) + 5-(r)(7+(r) 



rl3 
Jo 



(33) 



where E^j describes the self-consistently determined couphng with the conduction 
electron bath. Here cj.(r) and c^(t) are Grassmann variables fl^. Now DMFT for 
the DE model is exactly solvable in the classical spin limit 5* ^ oo where S(r) becomes 
r-independent (we scale J and h as 1/S), since for S constant 5* is diagonal in the 
Matsubara frequency representation: 

E„t + J/2 J/2 S- 1 / c„T ' 

J/2 S+ Eni - J/2 \ \ Cni 



S 



E 



C, 



nl 



(34a) 
(346) 



defining the matrix An and the components S'„ of S. Here Ena = — Jq dr exp(zco'„r)i?o-(r) 
and Un is a fermionic Matsubara frequency. Note that as shown in the self- 
consistency condition for Ena can be written as Ena = T,{iuJn) + G{iuJn)~^, so Ena- 
is just our quantity Erj{iujn). 
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S= J= oo , n= 0.75 



3.0 I 1 1 1 1 1 r 




0.0 0.2 0.4 0.6 0.8 

Temperature, kgl/W 



Figure 3. The S = J = oo inverse magnetic susceptibility x ^ for n = 0.75 compared 
witli tlie Curie law and DMFT values. 



In terms of S the partition function Z is given by 
Z-f d'S I (ndcLdc..)e-^ = / d-Se^^^' n / (ndcLdc..)e-^" 



(35) 



where / d^S is the integral over the local spin direction. All site-diagonal correlation 
functions can be calculated explicitly in terms of the EnaS, for example the one electron 
Green function is given by 



1 



(36) 



It is convenient to work with the generating function 



J (jl dcLdc^^j CXp (^-Sn + Yl (jia^na + cL^na)^ (37a) 



det(A„) exp 



( ) 



-1 / Vn] 



(376) 



In terms of Zn the partition function Z = / d^S* exp{f3hS^) (Hn ^n),,/s=o 
/ d^iS cxp(f3hS^) Hn det{An) and the local spin probability distribution function P(S) 
Z~^exp{phS') {UnZn}r,'s=o- Notc that {U{S)) = Jd'^S P{S)U{S) for any U. 
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Explicitly, our full Green functions are given in = cxd DMFT by 

d2 



1 



r]'s=0 



7?'s=0 



and these expressions are easily evaluated to give 



J 



E. 



J/2S' 



(E„t + J/2 S^) {Eni - J/2 S^) - J^/A (l - {S^ f) 



1 - {S'^y 



(38 a) 
(385) 

(39 a) 
(395) 



(E„t + J/2 S^) {Eni - J/2 S^) - J74 (1 - {S^f) I ■ 

Summing the S = 00 CPA expressions ( |20a| ) and ( p06| ) over a we obtain the analytic 
continuations (from the Matsubara frequencies) of ( |39aD and ( p9^ ) respectively, i.e. 
in the classical spin limit our CPA agrees with DMFT for Sa{\) and Tct(A). This is 
important since DMFT is known to be exact for dimension D = 00 and is perhaps 
the most natural local approximation for D finite. 

Similarly we can calculate the a = ± components of these Green functions within 
DMFT as 



(40 a) 
(405) 



. , ^^AEni-J/2S-^){E^,^ + J/2S-^)-J'/A[l-{S-^f 



n' 

J 1 



det(A„) det(A 
V ^ ^ ^ det(A„)det(A„0, 



" n' 

and = S-^—S^ etc. Now for a function g{z) analytic off the real axis fi'(^'-^ 
I[g], so applying this to the n' summations in (|40q|) we find 

{Eni - JS^/2){h^{S') + JS'h{S^)/2) - JV4(l - {S^f)h{S^ 



where 



(E„^ + J/2 S^) {E^^ - J/2 S^) - JV4 1 - (5 



E-^ 

{E^ + J/2 S^) {E^ - J/2 S^) - JV4 (1 - {Si' 

1 



(41) 



(42 a) 



(425) 



(E-f + J/2 S') {E^ - J/2 S') - JV4 [1 - {S' 

It may be shown using the sum rules (|17a|) and (|17^ ) that ((S'^)™(/i|(5'^) + 
JS'l2{S')/2)) = ((^")™n|) and J/2 ((S^)™(1 - {S'f)l2{S')) = -{{S')"'S-a+) for any 
m, so in we can replace /it(^^) + JS'l2iS')/2 with and J/2 (1 - {S'f)l2iS^) 
with —S~a~^. A similar procedure can be carried out for ( [40 5|) , and we obtain 

{Eni - J/2S') + J/2 S-a+ ^ 



5*1 (A, iujr, 



{Er,^ + J/2 S^) (Ki - J/2 S^) - JV4 (1 - (5 



(43 a) 
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/ X <^/2 (iS') - lUt - (Era - <^/2 S') S-a+ 

T+(A, ZLOn) = ( e''' ^ ^ ^ ^ 1 ^ ) . (43 &) 

^ \ (Kt + - J/2 5-) - JV4 (l - (S-^'' ^ 



It may be seen that the DMFT and CPA expressions for agree, but T|^(A)|dmft = 
T|"(A)|cpA- This discrepancy, due to failure of the decouphng approximations made in 
the equation of motion for ||^, vanishes in the important hmit J — *• cxo of strong 
couphng where all a = + Green functions are zero. 

5. The Curie temperature 



Furukawa [|T4| finds that the S = J = D = oo DE model exhibits a transition to 
ferromagnetism, and in section ^ it was shown that our CPA gives exact expressions 
for the Green function equations in this limit. However, in section ^ we showed that 
there is no magnetic transition in our CPA (figure |I]) and found unphysical behaviour 
at n = 1 (figure |^). In this section we resolve this apparent discrepancy and postulate 
a modification of our CPA that restores magnetic behaviour. We again work in the 
elliptic DOS case where En^- = iojn — W"^ / AG„{iuJn)- 

We first derive the S = J = oo DMFT static susceptibility x- In this case 

p ^ exp i(3hS^) Un [(1 + S^) E„T + (1 - S^) i^nj 

^ ' Jd^Sexp{(3hS^)nn[{^ + S^)E^^ + {l-S^)E^^y ^ ' 

as was shown in the previous section. The first order deviation 5-P(S) of -P(S) from its 
zero-field paramagnetic state value (l/47r) is given by 

where En is the zero-field paramagnetic state value and 5En is the first order deviation 
of En^. From (^) we have in the present case, where W = W/V^, 

^Gn. = 2E^^-wy 12- ^^^^ 

Substituting (^ into (^ using 6Ena = —W^/A5Gna, then multiplying by S"" and 
integrating over the local spin orientation we obtain 

9f^B5{S^) ^ {9^^BfP/^ ,4^. 

B 1 + W^V12 En(2^^- W^V12)"' 

upon rearrangement. This gives the correct Curie law x = 1/3 (s'/^b)^/^ at n = 
and n = 1 and a transition to ferromagnetism for all n at a temperature k-^T = 
-Wyi2I[{2E^ - Wyi2)-^] where E{e) = e - W^/AG{e) H. 

We now consider where the CPA calculation of x gone wrong and how to 
improve it. The CPA equations for the Green functions are exact in the present case 
{S = J = 00), but to solve them we need expressions for the expectations (exp(AS'^)no-), 
(exp(AS'^)S'^°"cr°"), and (exp(AS'^)), as mentioned in section ||. Since the first two of these 
are obtained using the I sum rule, a procedure that is exact, the problem must lie with 
the determination of (exp(AS'^)). Note that the way that we have calculated (exp(AS'^)) 
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only works for S finite, and to calculate it at 5* = oo we have worked for finite S and 
taken the limit at the end. 

Now knowledge of (exp(AS'^)) is equivalent to knowledge of -P(S), so a possible 
way of improving the CPA expression for x is to abandon the above self-consistent 
determination of (exp(AS'^)) and instead to use some expression for P(S) that reduces 
to the S = oo DMFT result (^4|) in the classical spin limit. We have so far been unable 
to derive such an expression, so we instead postulate a natural extrapolation of (^) to 
finite S, justifying our formula by the resulting behaviour of x- This procedure will at 
least force the CPA for x to become exact in the S oo limit, and the S = oo magnetic 
transition is likely to persist to finite 5*. 

From the S" ^ oo limit of (|T9^ and (|T9^ ) it may be seen that the quantity in 
square brackets in (|^) is related to the denominators of the Green function equations 
(for S = oo). A natural extension of (|4^ ) is thus 



p .g^ ^ exp iPhr^Sn Un [(1/2 + 3 + 8^ E^^ + (1/2 + S - S^) E^^] 

^ ' 27rE5^exp(/5/ir/5-)nn[(l/2 + 5 + 5-)i?„T + (V2 + ^-5^)Kj ^ ' 

where t] is chosen so as to optimise the behaviour of x- The explicit energy shifts 
associated with the field h in E-^ and E^ are ambiguous and have been neglected. 
However some effect of h on the double exchange enters through the Green functions and 
the factor t] allows for the conduction electron contribution to the spin in the Zeeman 
energy. Then proceeding as above we can calculate 

2 1/3 vf3S{S + 1) - (3sI[G'/E] + 1/2 1[G'] 
^ = ^^^^^ l-Q/i2S+l) + 2PsR ^^'^ 

where s = W^U S{S + 1)/{2S + 1), Q is as in (|9^, and 

1 



R 



(50) 



(25 + 1)^2 _ (2/3 S + l)WyA_ 

It is then easy to see that if we take r] = 1 + {n/2)/{S + 1) then the correct Curie laws 
X = {gnBypS{S + l)/3 and x = (gf^EfPiS + 1/2){S + 3/2) /3 are obtained at n = and 
n = 1 respectively. Note that St] = S + n/2 + 0{l/S), so if we regard our extrapolation 
from S" = oo as a kind of 1/S expansion this is a very natural value — it corresponds to 
the average spin size in the system to leading order. With this form ( ^^D for x we also 
obtain a magnetic transition for all 5* at a temperature determined by 

2sR , , 

The Curie temperature is plotted against filling n for various S in figure |^ (left figure) 
below. It agrees with Furukawa's result in his case of 5* = oo. Clearly for finite S 
ferromagnetism is more stable for n > \j2 than for n < 1/2, in agreement with the 



findings of Brunton and Edwards [T^. We have also calculated Tc via the spin- wave 
dispersion in the (assumed) saturated ferromagnetic groundstate, using a method similar 
to that of Sakurai [0, ^ for the Hubbard model. The Curie temperature obtained is 
similar in magnitude to Tc in figure ^ and decreases with increasing S", as is the case 
here for n near 1. This work will be published elsewhere. 
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Brunton and Edwards found that the stabihty of the spin-saturated state at T = is 
strongly dependent on the bare DOS used: approximating the true cubic tight-binding 
DOS with the elhptic DOS quahtatively changed the form of their spin-flip excitation 
gap. Accordingly we check the effect on Tc of using the true tight-binding DOS. The bare 
elliptic and cubic tight-binding DOSs and the corresponding full (zero-field paramagnetic 
state, S" = 1, n = 1/2, and J = oo) CPA DOSs are shown for comparison in figure ^ 
below. Now it is straightforward to extend the derivation of the Tq equation to the case 
of a general DOS; the only effect on (§T|) is to replace W^/A with {W/{WG)y + l/G' 
inside the I-functionals. Hence for general DOS 

2S(5+1) R 
3' 



where 



R 



Q 



25 + 1 g/(25+l)-l 

{W/(WG)f + l/G" 



(52) 



(25 + 1)^2 - (2/3 S + l)W^/W^{{W/{WG)y + l/G') 



E 



(53a) 



(536) 



- vm^/W^{{W/{WG)y + l/G') 

Note that for the elliptic DOS case {W/{WG))'^ + l/G' = W^A. We plot Tc obtained 
from (|5^) in the most sensitive case of 5 = 1/2 in figure ^ (right figure) for the elliptic 
and cubic DOSs. It may be seen that changing the form of the bare DOS does not 
have a large effect on Tc. The dip in Tc near n = 0.3 for the cubic DOS is interesting 



since near this filling Brunton and Edwards |T^ found an instability of the saturated 
ferromagnetic state for S = 1/2. 




-0.4 0.0 0.4 
Energy 



Figure 4. Bare (left figure) and full CPA (right figure, ioi S = 1, n = 1/2, J = oo, 
W = 1 zero-field paramagnetic state) DOSs. 



6. Summary and outlook 



In this paper we have extended our many-body CPA treatment of the DE model 
to the case of general 5* and magnetisation. In our original approach we were faced 
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Filling, n 



0.04 




0.00 1 • 1 ■ 1 • 1 • 1 • 1 

0.0 0.2 0.4 0.6 0.8 1.0 

Filling, n 



Figure 5. The Curie temperature k-QTc/W calculated using the elliptic bare DOS 
plotted against filling n for various S (left figure), and the effect on Tq for 5* = 1/2 of 
changing the bare DOS to the 3d-cubic DOS (right figure). 



with 45 + 1 algebraic equations to solve for the Green functions in the case of non- 
zero magnetisation. A correspondingly large number of correlation functions had to be 
determined self-consistently. Consequently in we only considered S = 1/2 for the 
magnetised state and subsequently [1^ calculated the paramagnetic susceptibility in 
this case. The generalisation to arbitrary S in section |^ of this paper is achieved by 
introducing generating Green functions involving a parameter A. The 45* + 1 coupled 
algebraic equations are then replaced by a single first order linear differential equation 
in A whose solution yields the CPA equations for the Green functions. Only three 
correlation functions have to be determined, as functions of A, and two of these may 
be obtained directly from the Green functions. The indirect determination of the third 
(exp(AS'^)), from the approximate EOM for the Green functions, is less reliable. It 
seems to be the origin of difficulties in section H, where the paramagnetic susceptibility 
is calculated for J = oo. No ferromagnetic transition is found for any n or S and for 
n = 1 the correct Curie law, with spin S + 1/2, is found only at high temperature. On 
the other hand in section ^ it is shown that for S = oo, where dynamical mean field 
theory has been implemented |14], our CPA equations for the Green functions agree 
with DMFT. Furthermore DMFT leads to a ferromagnetic transition for < n < 1 and 
to a correct Curie law for n = 1. In section ^ this paradox is resolved by abandoning 
the apparently unreliable self-consistent determination of (exp(A5'^)) and using instead 
a probability distribution P{S^) to evaluate the required expectation values. The form 
of P{S^) used for finite S* is a reasonable extension of the form which arises in DMFT 
for S = oo. We then find a finite Curie temperature Tq for < n < 1, and correct Curie 
laws for n = and 1, for all S. Naturally the results agree with DMFT for S = oo. 
The maximum in Tq, as a function of band-filling n, moves from n = 0.5 for 5 = oo to 
larger values of n as decreases. 

This work completes our present study of the paramagnetic state and ferromagnetic 
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transition of the DE model within our many-body CPA. With some effort we could 
pursue the calculations into the ferromagnetic state. However this has already been 
done for 5 = oo within DMFT and the rewards might be slight, particularly since 



for finite S the CPA never gives a ground state of complete spin alignment. It seems 
more profitable to repair some defects of the DE model itself. One should include both 
coupling to phonons and the double degeneracy of the eg band. It is likely, as originally 
proposed by Millis et al that phonon coupling is essential for an understanding of 
the insulator-like paramagnetic state in the manganites. We showed that, without 
phonons, the DE model gives much too small a resistivity. The introduction of phonons 
is therefore a high priority and it is in fact easier to include coupling to local phonons 
in our CPA approach than to consider degenerate orbitals. This is our next objective. 
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